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BAND-LIMITED ADAPTIVE FEEDBACK 
CANCELLER FOR HEARING AIDS 


BACKGROUND OF THE INVENTION 

1. Related Application 

This application claims the benefit of co-pending provisional application 
Serial No. 60/102,557 filed September 30, 1998. 


2. FIELD OF THE INVENTION 

The present invention relates generally to the field of audio amplification, 
especially for hearing aids. More particularly, the invention provides a method 
for efficiently cancelling acoustic feedback in the hearing aid. 

2. PRIOR ART 

Acoustic feedback in a hearing aid^ front the electroacoustic transducer 
(commonly referred to as the "receiver") back to the microphone, is common and 
difficult to suppress. Feedback may produce an audible whistle, which is 
irritating to the hearing aid wearer so that the wearer must often reduce the volume 
to a lower than desired level, thereby reducing the effectiveness of the hearing aid. 

One reason that it is difficult to effectively suppress acoustic feedback in a 
hearing aid is that the frequency at which feedback occurs varies with changing 
external conditions. Therefore, for feedback to be effectively cancelled without 
undesirably degrading the amplified signal, some form of adaptive cancellation is 
required. Various techniques have been proposed for implementing adaptive 


> 
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feedback cancellation. Such techniques are disclosed, for example, in U.S. Patent 
Nos. 5,016,280; 5,091,952 and 5,259,033. 

One of the principal design objectives for hearing aids is miniaturization of 
physical volume. Most wearers prefer a hearing aid that can be worn entirely 
within the ear. Advances in microelectronics have vastly improved the signal 
processing capabilities of in-the-ear (ITE) hearing aids. Even so, the provision of 
effective feedback cancellation remains a practical design challenge. Prior art 
techniques necessarily require certain trade-offs. As a result of such trade-offs, 
the hearing aid may exhibit only a small increase in maximum stable gain, slow 
filter adaptation, distortion, interference and/or lack of adaptation to individual 
wearers. 
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SUMMARY OF THE INVENTION 

The present invention provides an improved method for adaptively 
cancelling acoustic feedback in hearing aids and other audio amplification devices. 
Feedback cancellation is limited to a frequency band that encompasses all unstable 
frequencies. By limiting feedback cancellation in this manner, a relatively simple 
signal processing algorithm may be used to produce highly effective results with 
minimal signal distortion. 

In order to implement the present invention, unstable feedback frequencies 
must first be identified. This is accomplished by various techniques with real ear 
measurements, from which the complex open loop transfer function may be 
derived. Once the unstable feedback frequencies are identified, a band limited 
adaptive filter is implemented. By thus limiting the bandwidth of adaptation, the 
adaptive feedback canceller is able to adapt very quickly within the range of 
unstable frequencies with relatively low adaptation noise. By limiting the 
bandwidth of the feedback cancellation signal, the distortion due to the adaptive 
filter is minimized and limited only to the unstable feedback regions. Compared to 
broadband feedback cancellation, the band limited feedback canceller of the present 
invention produces less distortion, and therefore the output sound quality is much 
improved. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a functional block diagram of a hearing aid in which the present 
invention may be implemented. 

Figure 2 is a simplified diagram illustrating the transfer functions of the 
hearing aid of Figure 1. 

Figure 3 illustrates a prior art technique for deriving the open loop 
transfer function of a hearing aid. 

Figure 4 is a functional block diagram of a hearing aid which incorporates 
a generic broadband adaptive feedback canceller. 

Figure 5 is a functional block diagram of a hearing aid having an internal 
noise generator forderivation of the open loop transfer function. 

Figure 6 is a functional block diagram of one embodiment of a hearing aid 
having a band-limited adaptive feedback canceller in accordance with the present 
invention. 

Figure 71s a functional block diagram of another embodiment of a hearing 
aid having a band-limited adaptive feedback canceller in accordance with the 
present invention. 

Figure 8 illustrates base-2 logarithm quantization as used by the present 
invention for filter coefficient adaptation. 

Figure 9 is a functional block diagram illustrating the signal processing 
used to implement an adaptive digital filter. 

Figure 10 illustrates the processing used for power estimation. 
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Figure 11 illustrates the processing used for coefficient adaptation. 

Figure 12 illustrates the processing used for DC removal. 

Figure 13 illustrates the processing used for bandpass filtering. 

Figure 14 is a functional flow diagram illustrating the process scheduling 
implemented by the control unit of Figure 9. 

Figure 15 is a logic diagram corresponding to the flow diagram of Figure 

14. 

Figure 16 is a functional block diagram of an alternative embodiment of 
the present invention. 

Figure 17 is a functional block diagram of another alternative embodiment 
of the present invention. 

Figure 18 is a functional block diagram of still another alternative 
embodiment of the present invention. 

Figure 19 is a functional block diagram of yet another alternative 

embodiment of the present invention. 

Figure 20 is a functional block diagram of adjustable FIR filtering with a 
power-of-two scaling gain. 

Figure 21 is an example of open loop transfer function and unstable 
feedback frequencies measured on KEMAR ear. 
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DETAILED DESCRIPTION OF THE INVENTION 

In the following description, for purposes of explanation and not 
limitation, specific details are set forth in order to provide a thorough 
understanding of the present invention. However, it will be apparent to one 
skilled in the art that the present invention may be practiced in other embodiments 
that depart from these specific details. In other instances, detailed descriptions of 
well-known methods and devices are omitted so as to not obscure the description 
of the present invention with unnecessary detail. 

Figure 1 is a functional block diagram of a hearing aid 10 in an ear of a 
•hearing aid user. The hearing aid 10 comprises a microphone 12, microphone pre- 
amplifier circuitry 14, analog-to-digital converter 16, signal processing circuitry 
18, digital-to-analog converter 20 and receiver 22. The signal processing circuitry 
18 has a transfer function from point C to point D of K(f). The feedback path 
from point D to point C includes digital-to-analog converter 20, hearing aid 
receiver 22, acoustic and mechanical coupling between the receiver and 
microphone, analog conditioning circuitry 14, and analog-to-digital converter 16. 
P(/) is the transfer function of the feedback path. 

Figure 2 shows a simplified hearing aid model in the ear. The open loop 
transfer function of the hearing aid is defined as K(f)p(/). A hearing aid may 

produce unstable feedback only at those frequencies,/, , where its open loop 
transfer function meet the Nyquist instability criteria: 
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Magnitude: \K{f u ) P(f tt )\ > 1 

Phase: ZK(f u ){3(f u ) = nx360» 

Equation J 

When the above phase condition is met and the magnitude of the hearing aid open 
loop transfer function is close to but less than unit gain, the hearing aid is in a sub- 
oscillatory state. Both sub-oscillation and oscillation are unpleasant to hearing aid. 
wearers and must be eliminated. In the following description, We refer both 
unstable feedback frequency and sub-oscillatory frequency as unstable feedback 
frequency. An efficient feedback cancellation algorithm needs to operate only in 
the frequency region containing all f u . We have found that when feedback 
cancellation is limited to this frequency region, a relatively simple algorithm can 
produce highly effective results with minimal signal distortion. None of the 
previously developed methods of feedback cancellation in hearing aid systems has 
exploited this discovery. 

To derive the open loop transfer function^ conventional prior art methods 

break the closed loop system after the microphone pre-amplifier and before the 
hearing aid processing module, as shown between points C and D in Figure 3, and 
make a standard two-channel transfer function measurement with a signal 
analyzer. However; these methods cannot be conveniently used with commercial 
hearing aids, because it is impractical to break the hearing aid circuit. We have 
developed the following three methods to derive the open loop transfer function 
without breaking the hearing aid circuit. 
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The first method utilizes a probe microphone to accurately derive the open 
loop transfer function K(j)$(J) from a closed loop measurement by acquiring the 

probe microphone signals near the eardrum with the hearing aid in the ear canal. 
This method is described in co-owned patent application Serial No. 08/926,320, 
the specification of which is incorporated herein by reference. Briefly, the method 
requires acquiring probe microphone signals near the eardrum to compute the open 
loop transfer function with the hearing aid in three gain settings, Gl , G2, and G3, 
where G3=0 (i.e., the hearing aid is off). The method and procedure of detecting 
and digitizing the acoustic signal with a probe microphone used in the current 
invention are, in part, the subjects of co-owned U.S. Patent No. 5,325,436, 
"Method of Signal Processing for Maintaining Directional Hearing with Hearing 
Aids". 

The complex open loop transfer function AT(/)p(/)can be derived as: 

K(f)B(f) = Gl X //2ab( ^ ) Z G1 x H iab( /> ± Hlf^Jl x < G2 ~ G1 > 

G2 x Gl x (//2ab(/) - H fab(/)] 

Equation 2 

where H iAB (f), H 2A b(/) and H 3AQ (f) are closed loop transfer functions from point A 
to point B in Figure 1 measured with Gl, G2 and G3. Point A is at the hearing 
aid microphone input and point B is at the hearing aid receiver output, i.e., at the 
probe microphone. 

Because the probe microphone must be used for the closed loop 
measurements, the venting system of the hearing aid is partially blocked by the 
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probe tube. The derived open loop transfer function is thus only an 
approximation of the open loop transfer function with an open hearing aid vent. 
The closed loop transfer function measurements are also sensitive to head 
movement and surrounding environmental changes during the measurements. 

The second method does not require a probe microphone measurement. It 
requires configuring the adaptive feedback canceller in the hearing aid to operate in 
a broadband mode as shown in Figure 4. The adaptive digital filter (ADF) 30 in 
the hearing aid provides an estimation of the feedback path including the digital-to- 
analog converter 20, the receiver 22, the coupling path between the receiver and 
microphone, the microphone 12, the AGC 14, and the analog-to-digital converter 
16. 

Data is collected with the hearing aid wearer seated in a quiet room. A 
white noise signal is generated in the sound field through a loudspeaker. The 
hearing aid is programmed with a known reference response. The hearing aid 
wearer is instructed to adjust the hearing aid gain below the Uncomfortable Level 
(UCL). While the white noise is playing in the sound field, the ADF in the 
hearing aid adapts itself to match the feedback path. After the ADF has 
converged, the filter coefficient vector W(w) is read from the hearing aid together 
with the "delay". The feedback transfer function p'(f) can be estimated as 
follows: 
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P'(f) = e-^xW n (f) 

Equation 3— 

where D is the delay introduced by "Delay" block 32; W n {J) is the frequency 
response of the ADF at time index n, for which the impulse response is W(«). 
To improve the accuracy of the measurement, several coefficient vectors W(w,) 
with different time index n f , can be obtained from the ADF and averaged to 
compute W n (J) in Equation 3. 

In order to compensate for the hearing loss of each individual, the hearing 
aid reference response used during the measurement is replaced by a desired 
hearing aid response K(J). The complex open loop transfer function associated 
with the desired hearing aid response K{f) can be estimated by K(f)p'(f)\ 

K(f)P , (f) = e- j27rJD xW n (f)xK(f) 

Equation 4 

The unstable feedback frequencies can then be determined by the estimated 
open loop transfer function based on Equation 1. Comparing with the first 
method, this method does not require a probe microphone measurement. 
Therefore, the venting system of the hearing aid is not blocked during the 
measurement. The measurement is not sensitive to head movement and 
surrounding environmental changes. It does not require additional circuitry in the 
hearing aid to support the measurement. The white noise signal serves as a test 
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signal for the measurement,. but it also acts as interference noise to the ADF since 
the signal directly from the loudspeaker does not carry any information about the 
feedback path. Thus, it introduces adaptation noise to the ADF coefficients. To 
achieve the best results, the level of the white noise signal should be kept as low 
as possible but above the noise floor of the room. The hearing aid gain should be 
set as high as possible but below UCL and without audible feedback. 
Furthermore, the hearing aid is configured in a closed-loop during the 
measurement. The unstable feedback rather than the UCL may be the limiting 
factor for setting the hearing aid gain. 

The third method is similar to the second method discussed above. An 
internal pseudo-random noise generator as shown in Figure 5 is used to replace 
the external white noise test signal, and the hearing aid processing block is 
disconnected. This method has two advantages over the second method discussed 
above: (1) the adaptation interference introduced by the external test signal is 
removed by using an interna! pseudo-random noise generator; and (2) the hearing, 
aid is operated in an open-loop configuration during the measurement. Thus, the 
hearing aid gain can be set as high as possible without instability problems. This 
method requires adding a noise generator and some control logic to the hearing aid. 
The open loop transfer function measurement is only required during the initial 
fitting process. 

Figures 4 and 5 both illustrate a generic broadband adaptive feedback 
canceller. A digital feedback cancellation signal is generated by passing the hearing 


W ° 00/19605 PCT/US99/22757 

12 

aid digital output signal at point B to an ADF 30, which approximates the 
feedback path p(#. The digital feedback cancellation signal is subtracted from the 

digitized microphone output at point C, eliminating the feedback signal from the 
hearing aid. 

The ADF 30 comprises an adjustable filter, which uses filter coefficients 
to generate a feedback cancellation signal. A coefficient adaptation portion 34 
adjusts the filter coefficients to approximate the feedback path. There are various 
adaptive filtering methods, which have different filtering structures and different 
adaptation algorithms. Some of them exhibit extremely good convergence behavior 
and accuracy, but require extensive computation. For example, the family of 
Least-Squares (LS) algorithms belongs to this category. For the application of 
adaptive feedback cancellation to hearing aids, the simple Stochastic-Gradient 
(SG) algorithms are sufficient to provide acceptable performance. The simplicity 
of the adaptive filtering algorithm is a very important factor for hearing aid 
applications since it is desirable to minimize the hardware requirements. We 
choose an FIR (Finite-Impulse-Response) filter with Normalized Least-Mean- 
Square (LMS) algorithm for the ADF coefficient adaptation because of its 
simplicity and satisfactory convergence performance. To make the adaptive 
filtering algorithm more specifically suitable for the problem of adaptive feedback 
cancellation in hearing aids, we adopted a Normalized LMS algorithm with the 
adaptation step-size normalized to ADF input and error, as discussed below. 
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In general, the adaptive filter is a time-varying system, which can track the 
time- varying feedback paths. But unnecessary variations due to the adaptation 
noise of the adaptive filter may also introduce unpleasant distortion in the hearing 
aid output. Thus, the adjustment of the adaptive filter should be minimized except 
as necessary for eliminating the unstable feedback. 

The signals coming from external sound sources (external input signals), 
other than the signals which are fed back from the hearing aid output to the hearing 
aid microphone, are usually considered as an interference to the adaptation of the 
adaptive filter. An external input signal normally introduces bias to the optimal 
solution for the adaptive filter. The bias varies when the external input signal's 
property changes. The external input signal also causes misadjustment noise 
during the adaptation process. By removing the frequency components that are 
unlikely to have potential oscillation problems (unstable and sub-oscillatory) from 
the adaptation signals: e(n) and x(n) 9 which are used for the adaptation process , 
the interference effects of the external input signal on the adaptation process can 
be greatly reduced. 

For the modified Normalized LMS algorithm, the convergence speed of the 
adaptive filter is governed by the adaptation step-size, which is inversely 
proportional to the combined signal power of e(n) and x{ri). Reducing the 
combined signal power would lead to a larger adaptation step-size, which in return 
increases the convergence speed. Since we are removing the components that are 
irrelevant to the oscillation problems from the adaptation signals, the combined 
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signal power of e(n) and x{n) is reduced so that a larger adaptation step-size may 
be used to increase the convergence speed. With a higher convergence speed, the 
feedback canceller can better track the dynamic feedback path and the sudden gain 
change of AGC due to a sharp transition of input signal. 

Furthermore, due to the recursive topology used in the adaptive feedback 
cancellation, the periodicity in the external input signal would cause the 
cancellation of original signal (external input signal)at the hearing aid output. By 
removing the frequency components that are irrelevant to the oscillation problems 
from the adaptation signals, the periodicity in the adaptation signals, which 
originates from the external input signal, can be reduced. Therefore, the problem of 
canceling the original input signal at the hearing aid output will be relieved. 

In addition, the adaptive filter always functions better at frequencies where 
large energy exists. This means that the adaptive feedback canceller actively works 
at oscillation frequencies only when the energy of the oscillation components 
existing in the adaptive filter input signal and the error signal is comparable or 
greater than the peak energy of the spectrum of external input signal. As a result, 
the magnitudes of the frequency components at unstable feedback frequencies are 
build up and then suppressed, up and down, around the levels of the external 
input spectral peaks. This is so-called modulation effect of residual oscillation 
components. 
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By removing the unnecessary frequency components from the adaptation 
signals, the peak spectral level is much reduced. Therefore, the magnitudes of the 
residual oscillation components are also significantly reduced. 

This concept of removing unnecessary frequency components from the 
adaptation signals is extremely important and beneficial when the external input 
signal is a speech signal. As is well known, speech contains most of its energy 
and periodicity at the low frequencies, at which unstable feedback is unlikely to 
occur. Therefore, the low frequency components of speech are considered 
unnecessary for the adaptive filtering and can be easily removed by applying a 
highpass (or bandpass) filter to the adaptation signals. The cutoff frequency of the 
highpass filter is normally set at 200 Hz below the lowest unstable feedback 
frequency. 

As indicated above, the unstable feedback only occurs at certain 
frequencies. To effectively suppress oscillation at these frequencies without 
~ introducing distortion at other frequencies, the current invention configures the 
adaptive feedback canceller in such a way that it limits the bandwidth of 
adaptation signals to the frequency regions known to contain the oscillation 
frequencies. By doing so, the adaptive feedback canceller adapts very quickly in 
the oscillation frequency regions with much less adaptation noise and adapts very 
slowly in other regions. As a result the feedback cancellation signal generated by 
the ADF is also limited to the same frequency regions. Therefore, unlike the 
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broadband scheme, the band-limited feedback canceller produces less distortion, 
and hence the output sound quality is much improved. 

Figure 6 shows one possible structure to implement the band-limited 
adaptive feedback canceller. BPF1 and BPF2 are bandpass filters, which only 
pass the frequency region containing all the unstable feedback frequencies. BPF1 
and BPF2 remove most of the unnecessary frequency components from e(n) and 
x(n) to improve the adaptation process. The hearing aid input signal is filtered by 
BPF2 and then used as the desired signal for the adaptive filtering so that the 
adaptive filter only generates the cancellation signal in the passband of BPF2. A 
A 2 -delay linear-phase FIR filter is selected as the bandpass filter BPF2 so that 


we 


can use a pure delay A 3 to ensure the unfiltered hearing aid input d(n-A 2 ) is 

properly time-aligned with the filtered desired signal cf(n) at the cancellation 
junction SUM2 without introducing phase distortion in the primary hearing aid 
signal path. Since the additional delay A 2 in the primary hearing, aid signal path 
may increase the number of oscillation frequencies and shift the oscillation 
frequency to a lower frequency, the filter length of BPF2 must be minimized. The 
A, -sample delay in front of filter BPF1 is used to ensure the feedback cancellation 
signal is properly time-aligned with the feedback signal so that the impulse 
response of the adaptive filter contains most of the energy of the estimated 
impulse response of the feedback path. On the other hand, the delay introduced 
by the filter BPFI must not be too long or the cancellation signal generated by the 
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adaptive filter may lag behind the feedback signal. The group delays introduced 
by each block must meet the following condition: 

Da, (/)+ Z>BPFl(/) + DaDF_FR (f )= Df^dHurkpcth(f )+ £>BPF2 (f ) 

Equation 5 

The optimal delay Aj in samples can be determined based on the measurement of 

the feedback path during the hearing aid fitting process. In order to utilize all the 
ADF coefficients, the delay caused by the adaptive filter should be minimized. 
The delay A { must meet the following condition: 

Da,(/J< D/^dbackp tll h(/yhDB?F2{f)-DBPF\(f ) 

Equation 6 

We choose: 

Aj = min{D w ^(/)+D BPR (/)-0 Hm (/)} / ^^-e 

Equation 7 

Where e is the tolerance in samples. We typically choose e equal to 2. 

The output of the ADF is calculated one sample at time. At time index n y 
the calculation for the M-tap ADF output is described as: 
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M-l 


y(n)= y £w lt (n)-x(n-k) 


Equation 8 

where {w Q (n), w,(«), .... w Ml (n)} are the coefficients of the M-tap ADF at time 
index n. 

The coefficients of the ADF are updated with a modified Normalized LMS 
algorithm. As in the conventional Normalized LMS algorithm, the adaptation 
step-size is reduced when the signal level is high and vice versa. In feedback 
cancellation for hearing aids, however, the input signal of the ADF is delayed by 
the hearing aid processing module, BPF1 and A,. The power of the error signal 
e(n) is included for calculating the time-varying step-size as follows: 

e(n) = d'{n)- y ( n ) ^ where '(") is the output of the bandpass filter BPF2. 

Equation 9 

p(n) = p- p{n -l) + x 2 (n) + e 2 («) 


Equation JO 


M«) = 


[Q-p)-p(n) + ME]M 


Equation J J 


where p is the forgetting factor; 

C is a constant to control the adaptation speed; and 
ME is a small number to prevent the singularity of \x(n). 
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The coefficients of the ADF for the next time index w+l are updated as follows: 
w k (n+l)=w k (n) + ii(n)*e(n)'x(n-k) , for 0<k<M~l 

Equation 12 

Due to the bandpass characteristics of BPF1 and BPF2, the adaptation signals 
e{n) and x(n) only contain the frequency components in the oscillation frequency 
regions. 

Figure 7 shows another possible structure for implementing an adaptive 
band-limited feedback canceller. BPF1 is mainly used to limit the bandwidth of 
the cancellation signal and to limit cancellation artifacts to the cancellation 
bandwidth. BPF2 and BPF3 are used to limit the bandwidth of the adaptation 
signals to the' oscillation frequency regions. BPF1 , BPF2 and BPF3 are not 
necessarily implemented as linear-phase FIR filters, but their passbands must 
cover all the oscillation frequencies. 

The filtered adaptation signal samples are used for updating the time- 
varying step-size and ADF coefficients as follows: 

p(n) = p • p(n - 1) + x' 2 (n) + e' 2 f n; 

Equation 13 

_C 

A7-p;-pfn; + ME7-M 

Equation 14 
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w k rn + 7; = w k rn; + /iCn;eYn;xYn-k; , for 0<k<M-7 

Equation 15 

where x' (n) is the output of the bandpass filter BPF3; and 
e' (n) is the output of the bandpass filter BPF2. 

To make the adaptive system stable, the phase responses of BPF2 and 
BPF3 must be as close as possible. For simplicity and better adaptation, we select 
two identical high pass filters for BPF2 and BPF3. Unlike what is required for 
bandpass filter in Structure 1, the passband ripples of BPF1, BPF2, and BPF3 are 
not critical for good performance as long as the stopband attenuation is sufficient. 
In our experience, 30 dB stopband attenuation is adequate. Therefore, low-order 
IIR filters such as 2 nd or 3 rd order Elliptic IIR filters may be used for this 
application to reduce the hardware and computation complexity. 

In contrast to the embodiment illustrated in Figure 6, the various group 
delays in Figure 7 must meet the following condition: 

£>a, (J )+ DsrfAJ )+ Dadf_fir(/,)= D D Ac{f ;+fl K », ir W1<rf)Mt ( i /J + Dadc(J ' ) 

Equation 16 

The optimal delay A, in samples can be obtained in the measurement of the 

feedback path during the hearing aid fitting processing. In order to utilize all the 
ADF coefficients, the delay caused by the adaptive filter should be minimized. 
The delay A, most meet the following condition: 
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DAj {/ )< DDAc(f)+ Dacousiic feedback path{f)+ DaDc(/ )- DBPFl(/) 

Equation 17 

We choose: 

Aj = min{D DAC (/) + D^^Jf) + D ADC (/) - D.Jf)}^ " £ 

Equation 18 

Where e is the tolerance in samples. We typically choose e equal to 2. 

This band-limited feedback cancellation structure does not introduce any 
additional delay in the primary signal path and does not introduce additional phase 
distortion to the hearing aid output. 

The purpose of the ADF is to estimate the feedback path. In Figure 7, 
BPF1 is used to limit the bandwidth of the feedback cancellation signal. Because 
the frequency response of the feedback path has band-pass characteristics as 
generally shown in Figure 21, we may use BPFI to approximately match the 
frequency response of the feedback path. In this way, the ADF can be used 
mainly to track the variation of the feedback path. 

The band-limited adaptive feedback canceller can be implemented on a 
platform with either a general-purpose digital signal processor or a specialized 
digital signal processor. Due to the size and power constraints of hearing aid 
circuit design, it is desirable to utilize a fixed-point digital signal processor with 
limited precision and word length as the adaptive feedback canceller. Thus, the 
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efficient digital realization of the band-limited adaptive feedback canceller is 
extremely critical for the performance of the feedback cancellation under the 
constraint of limited hardware and computational resources. The present 
invention simplifies the computational requirements and addresses issues 
associated with the limited-precision effects on the adaptation process. 

In both of the above-described structures of band-limited feedback 
cancellers, the operations of adaptive filtering are performed in the same way. For 
a fixed-point implementation, some additional modules are required to maintain 
efficiency of calculation under the limited precision constraint. A generalized 
structure for adaptively filtered, band-limited feedback cancellers is shown in 
Figured. This generalized structure is applicable for both of the above-described 
embodiments. There are three input ports (*,(«), x 2 (n), and e(n)) and one output 
port (y(n)) for this adaptive filtering configuration. 

For the feedback canceller structure shown in Figure 6: 

x \ (") = x z(n) = x(n) 
For the feedback canceller structure shown m Figure 7: 

x,(n) = x(n) 
x 2 {n) = , x'(n) 
e(n) = e'(n) 

The generalized structure contains an adjustable FIR filtering module, a 
power estimation module, a coefficient adaptation module, a DC removin° 
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module, a coefficient bandpass filtering (CBF) module and a control unit. The 
purposes and detailed implementation of these modules is described below: 

1. Adjustable FIR filtering _ 

The adaptive FIR filter is used to approximate the dynamic feedback path 
and generate a feedback cancellation signal by convolving the input signal x { (n) 

with current filter coefficients {w k (n) ; 0<k< M-/}. However, the feedback 

path response has a very large dynamic range. In the fixed-point implementation, 
to fully make use of the internal precision of fixed-point adaptive FIR filtering, we 
must maximize the ADF coefficients to fit the word-precision allocated to them. 
An ADF scaling gain (G) is used to maximize the ADF coefficients and provide a 
wide dynamic range for the feedback cancellation signal. Therefore, the calculation 
of the A/-tap adaptive FIR filter output is slightly modified as: 
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k=0 


Equation 19 

The scaling gain G is selected as a power-of-two number 2 and can be 
implemented by left/right shifting. Normally, it is sufficient for L to be in the 
range [-3, 3], which provides a dynamic range from -1 8 dB to 1 8 dB. Figure 20 
shows a functional block diagram of adjustable FIR filtering with a power-of-two 
scaling gain. 

2. Coefficient Adaptation 

With the ADF scaling gain included, the calculation for time- varying step- 
size has to be modified as: 


Equation 20 

As shown in Equation 10 and Equation 13, p(n) is the power estimation of 
the combined signal e(n) andx 2 (/z). In the fixed-point implementation, because of 
the word-truncation in calculating pp{n-\ \ limit cycles will prevent p(n) from 

becoming zero, and therefore ME in Equation 20 is not necessary. For simplicity, 
we set ME to 0. 


WO 00/19605 


PCT/US99/22757 


25 


In order to avoid the division in calculation of time-varying step-size, we 
use a power-of-two number to approximate p(n) and select C, (1-p), G, and Mas 

power-of-two numbers: 


C = 2~ K 
(l~p) = 2- J 
G = 2 L 
M = 2 F 

Equation 20 can then be rewritten as: 


AT is a positive integer to control the adaptation speed. The range of K is 
typically from 7 to 10. The smaller value of K provides faster adaptation speed. 
J is a positive integer to control the time-constant of power estimations 
Typically, we choose 6 for J. L is an integer to control the ADF scaling gain. As 
stated above, the range is from -3 to +3. L is determined based on the feedback 
measurement during the hearing aid fitting process so that the filter coefficients of 
adaptive filter are maximized. F is an integer related to the length of the adaptive 
filter. With a sampling rate of 1 6000 Hz, we choose F = 5 (32 taps) so that the 
duration of the adaptive filter impulse response is about 2 ms, which is long 
enough to cover the variation of group delay in the feedback path in the unstable 
feedback frequency regions. ~ 



Equation 21 
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Q[ ] is a truncation operation. Q[log 2 (p(«))] can be implemented by 
searching the position index of the most significant bit (MSB) of pin). 

Figure 8 shows a functional block diagram of base-2 logarithm 
quantization. £ is a positive quantity represented in an unsigned binary integer 

format. The position index of the least significant bit (LSB) is 0. Q(log 2 (£)) 

returns the position index of the MSB of When £ is 0, Q[log 2 (£)] returns 0. 

We can further simplify the power estimation and coefficient adaptation 
by using a power-of-two number to approximate the error signal. This 
quantization of error signal does not affect the performance of the adaptive filter. 
Thus, the power estimation becomes: 

pin) = (1 - p) . p{n - 1) + x 2 2 {n) + e 2 (n) 

- p(n - 1) - 2~ J - p(n -1)4- x;(n) + sign(*(n)) • 2**** . 

Equation 22, 

The coefficient adaptation becomes: 

^(#1 + 1) = w k (n) + fi(n) ein) x 2 (n-k) , /or 0<A:< M-l 

- w, („) + 2^*-^"* W»l ■ sign(,(«» . • * 2 (n - it) 

= w^n) + sign(e(n)) ■ 2 y - |: - L - F - G f ,0 ^^ (n »M to ^^">»] - *) 

Equation 23 

Equation 22 shows that only one multiplication is required for the power 
estimation. The coefficient adaptation process shown in Equation 23 becomes a 


WO 00/1 9605 PCT/US99/22757 

27 

multiplication-free process and can be implemented with shifting, negation and 
addition operations. 

Figure 10 shows the functional block diagram for power estimation, and _ 
Figure 11 shows the functional block diagram for a multiplication-free, coefficient 
adaptation process in which: 

P(n) is the base-2 logarithm quantization of the error signal e(n); 

oc(n) is the base-2 logarithm quantization of the power estimation p(n); 

and 

v (n) is the sign of the error signal e(n). 

Since the band-limited feedback canceller works so efficiently and 
effectively in the unstable frequency regions, the adaptation step size used for the 
Normalized LMS adaptive filter may be further reduced to minimize the 
misadjustment noise and hence improve the hearing aid sound quality. With fixed- 
point implementation, the normalized LMS adaptation with a very small step size 
may be simplified to a sign LMS adaptation. Equation 23 may be rewritten as: 

w k (n + 1) = w t {n) + /usign(e(n) • x 2 {n - *)), for 0<k<M-l 

Equation 23a 

Where fi is a constant. For example, we may choose fi equal to 1 when w k (n) is 
represented with a 12-bit integer. By doing so, the computations associated with 
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power estimation, MSB search, and shifting that are required for the Normalized 
LMS are eliminated. 

3. Limited-precision effects due to fixed-point 
implementation 

In the fixed-point implementation of the adaptive-filtering algorithm, both 
the inputs and internal algorithmic quantities must be quantized to a certain 
limited precision. These quantization errors may accumulate without bound until 
overflow occurs, resulting in unacceptable performance. For example, a slight DC 
offset in e(n) and x 2 (n), which results from either the original ADC output or the 
word-truncation of band-limiting filtering, may accumulate over time and cause an 
increase of DC offset in the adaptive filter coefficients. The truncation operations 
in the adaptive filter coefficient adaptation may also cause similar DC build-up in 
the adaptive filter coefficients, especially when the signal level is low. 
Furthermore, the low frequency gain of the ADF filter response may gradually 
build up if the chosen band-limiting filter BPF1 has excessive low frequency 
attenuation compared to the feedback path. In both cases, the adaptive filter 
coefficients may overflow or saturate. A DC removing module is included to 
periodically remove the DC offset from the adaptive filter coefficient. Another 
bandpass filtering module is provided to filter the adaptive filter coefficients in 
order to suppress the low frequency and high frequency response build-up in the 
adaptive filter response. This operation is only needed when the filter coefficient 
is saturated. 
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In the DC removing module, the following operation is implemented to 
estimate the DC offset in the filter coefficient and subtract the estimated DC 
offset from the ADF filter coefficients: 


m(n) = — Y w k (n) 


Equation 24 

w k (n + 1) = w k (n) - m(n) , 0 < k < M - 1 

Equation 25 

Figure 12 illustrates the computational process for DC removal. 

For example, the DC removal may be scheduled every 256 samples at 
16000 Hz sampling rate. 

In the CBF module, a zero-delay bandpass filtering operation on the ADF 
coefficients is performed as follows: 

vv^(/? + l) = w^(w)/2-w Jt+l («)/2 > 0<A< M -3 

Equation 26 

w k {n+\)= w k (n)/2 , k = M - 2, M-\ 

Equation 27 

Figure 13 illustrates the computational process for bandpass filtering. 
This simple bandpass filtering operation doesn't offer a perfectly flat 0 dB 
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magnitude response across the passband, and may introduce minor audible 
distortion at the hearing aid output. It cannot be applied frequently and is only 
triggered when any one of the ADF coefficients is saturated. 

As mentioned earlier, the adaptive filter coefficient adaptation must be 
performed in company with DC removing and coefficient bandpass filtering 
operations. Since the DC removing and coefficient bandpass filtering operations 
don't need to be performed frequently, we schedule only one of these three 
operations to be performed at each sample period. Figure 14 is a flow chart of 
the scheduling process. In this way, the DC removing operation is performed 
periodically, and the coefficient bandpass filtering operation is triggered only 
when the SAT signal from the coefficient adaptation module is on. 

Figure 15 is a logical diagram of the control unit, which generates the 
control signals for the coefficient adaptation module, the DC removal module, and 
the coefficient bandpass filtering module. 

Experimental Results 

The first test of this invention was performed with a computer simulation. 
The simulation model was developed in SIMULINK and built with a dynamic 
feedback path. The dynamic feedback path was measured on a KEMAR ear with 
a clipboard slowly moving toward and away from the ear. Various hearing aid 
responses for human subjects were used as hearing aid processing for the 
simulation. 
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The hearing aid response and the dynamic feedback path were used to 
derive the open loop transfer function and identify the unstable feedback 
frequencies. The unstable feedback frequencies were used to configure the 
feedback canceller, specifically the bandwidth of the band-limiting filters. The 
tests were made with and without the band-limited adaptive feedback canceller. 
The maximum stable hearing aid gain was recorded under both conditions. The 
simulated hearing aid outputs were also used for subjective evaluation. 

The same tests were performed on human subjects with a real-time 
prototype digital hearing aid. The open loop transfer function of the hearing aid 
was determined based on closed loop probe tube measurement. The unstable 
feedback frequencies were identified from the open loop transfer function and 
used to configure the band-limited adaptive feedback canceller. The maximum 
stable insertion gains were recorded with and without the adaptive canceller. 

The results have shown that the band-limited feedback canceller effectively 
eliminates the oscillation and sub-oscillatory feedback and increases the stable 
hearing aid insertion gain by 12-15 dB with minimum distortion to the sound 
quality. 

Figure 21 shows an example of open loop transfer function measured on 
KEMAR ear. There are three unstable feedback frequencies as indicted by cross 
symbols. They are at 2660Hz, 3260Hz, and 3960Hz. Therefore, the cutoff 
frequency for the band-lirnitmg filters is set to 2460 Hz, which is about 200 Hz 
below the lowest unstable feedback frequency. 
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ALTERNATIVE EMBODIMENTS 

The embodiments shown in Figure 6 and Figure 7 are two specific 
examples of band-limited adaptive feedback cancellers. However, the invention 
embraces other embodiments having the same functionality in terms of limiting the 
bandwidth of adaptation and bandwidth of the cancellation signal to the frequency 
regions known to contain oscillation frequencies. 

For example, in Figure 6, we can replace the A r sample delay with two 
strictly complementary filters. One filter would be the same as BPF2, whose 
frequency transfer function is denoted as tf 2 (/),and the other would be its strictly 
complementary filter SCF2, which has a complementary frequency transfer 
function - H 2 (f). Such an embodiment is illustrated in Figure 16. 

We can further combine two identical BPF2s and split the hearing aid 
processing into two processes as shown in Figure 17. Hearing Aid Process 1 is 
used to process signals in the frequency regions covered by BPF2 T and Hearing 
Aid Process 2 is used to process the signals produced by BPF2's complementary 
filter, SCF2. Since the output bandwidth of hearing aid process 1 is the same as 
BPF2, the need for BPF1 as shown in Figure 6 is eliminated. A linear-phase FIR 
filter is selected for BPF2 and its strictly complementary filter is selected for 
SCF2 in order to minimize phase distortion in the primary signal path of the 
hearing aid. For some hearing aid applications, the phase distortion in the primary 
signal path can be either tolerated or corrected by other hearing aid processing 
modules. In this case, we can relax the strictly complementary condition to other 
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kinds of complementary conditions, e.g., power complementary so that infinite- 
impulse-response (HR) filters may be used to replace BPF2 and SCF2 to further 
reduce the hardware and computational complexity. 

So far, we have assumed that an analog automatic gain control (AGC) is 
used in the hearing aid between the hearing aid microphone and the analog-to- 
digital converter. Since the AGC is part of the feedback path but not part of the 
cancellation signal path, the sudden gain change of AGC due to a sharp transition 
of the hearing aid input signal may degrade the performance of the cancellation. 
One alternate construction, as shown in Figure 18, is to move the AGC after the 
feedback cancellation junction so that the AGC is in both the feedback path and 
the cancellation path. This would require a digital implementation of AGC. 
Another alternate solution, as shown in Figure 19, is to digitally apply the AGC 
gain to the feedback cancellation signal so that the feedback cancellation signal 
tracks the AGC gain change. This method would require converting the AGC gain 
to a digital format. 

The band-limited feedback canceller offers the best sound quality when the 
cancellation bandwidth is minimized. As is known, the oscillation frequencies of a 
hearing aid are functions of the hearing aid gain, which is normally controlled by 
the hearing aid volume control. Thus, the cancellation bandwidth of the feedback 
canceller can also be controlled by the hearing aid volume control in order to 
achieve the best sound quality. This can be accomplished by storing several sets 
of filter coefficients for the band-limiting filter in the hearing aid. With any given 
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hearing aid volume control setting, the appropriate set of filter coefficients is used 
to provide a filter response that covers all the oscillation frequencies and has the 
minimum bandwidth. This filter selection process is needed only when the 
hearing aid wearer changes the volume control setting. 

We may also choose different constant K to control the adaptation speed 
for different hearing aid volume control settings, since the higher hearing aid gain 
setting may require a faster adaptation speed to suppress the feedback. 

It will be recognized that the above-described invention may be embodied 
in other specific forms without departing from the spirit or essential 
characteristics of the disclosure. Thus, it is understood that the invention is not to 
be limited by the foregoing illustrative details, but rather is to be defined by the 
appended claims. 
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CLAIMS 

WHAT IS CLAIMED IS: 

1 . A feedback canceller for an audio amplification device comprising: 
an adaptive filter; 

means for combining an output of the adaptive filter with an input 
of the audio amplification device; 

a first band limiting filter having an input coupled to an output of 
the audio amplification device and an output coupled to an i&pttt^fjhe^r' 
adaptive filter, wherein the first band limiting filter has a passband limited 
to a frequency band containing unstable frequencies. 

2. The device of claim 1 wherein the audio amplification device 
comprises a hearing aid amplifier. 

3. The device of claim 1 wherein the adaptive filter is implemented in 
a digital signal processor. 

4. The device of claim 3 wherein the adaptive filter is a finite impulse 
response filter. 

5. The device of claim 4 wherein the adaptive filter is implemented 
with a normalized least mean square algorithm. 

6. The device of claim 1 wherein the first band limiting filter is 
implemented in a digital signal processor. 
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7. The device of claim 1 wherein the first band limiting filter 
comprises a high pass filter. 

8. The device of claim 7 wherein the passband of the high pass filter 
has a cutoff frequency approximately 200 Hz below a lowest of the unstable 
frequencies. 

9. The device of claim 1 wherein the adaptive filter comprises means 
for convolving the input of the adaptive filter with filter coefficients. 

10. The device of claim 9 wherein the adaptive filter further comprises 
means for applying a scaling gain to the filter coefficients. 

1 1 . The device of claim 10 wherein the scaling gain is a power of 2. 

12. The device of claim 1 1 wherein the scaling gain is in the range of 2~ 3 

to2 3 . 

13. The device of claim 9 further comprising: a second band limiting 
filter having an input coupled to the input of the audio amplification device and an 
output; a third band limiting filter having an input coupled to the output of the 
first band limiting filter and an output; 

wherein the second and third band limiting filters have passbands limited to the 
frequency band containing unstable frequencies and the filter coefficients are 
functions of the outputs of the second and third band limiting filters. 
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14. The device of claim 13 wherein the second and third band limiting 
filters have matching phase responses. 

15. The device of claim 13 wherein the second and third band limiting 
filters have substantially identical characteristics. 

1 6. The device of claim 1 further comprising means for delaying the 
output of the audio amplification device prior to the input to the first band 
limiting filter. 

17. The device of claim 16 wherein the output of the audio 
amplification device is delayed by an amount of time that is a function of a 
feedback path delay. 

18. The device of claim 1 further comprising an automatic gain control 
(AGC) coupled between the means for combining and the input of the audio 
amplification device. 

19. The device of claim 1 further comprising an automatic gain control 
(AGC) coupled to the input of the audio amplification device, said AGC having an 
output combined with the output of the adaptive filter. 

20. The device of claim 1 wherein the adaptive filter has an adjustable 
adaptation speed. 

21 . The device of claim 20 wherein the adaptation speed is adjusted as 
a function of a gain of the audio amplification device. 
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22. The device of claim 1 wherein a length of the adaptive filter is 
selectable. 

23. The device of claim 22 wherein the length of the adaptive filter is 
selected as a function of sampling rate and the frequency band containing unstable 
frequencies so that duration of the adaptive filter's impulse response is long 
enough to cover an expected range of variation of group delay in the frequency 
band containing unstable frequencies in a feedback path from the output of the 
audio amplification device to the input thereof. 

24. The device of claim 1 wherein the first band limiting filter has a 
frequency response approximately matching a frequency response of a feedback 
path from the output of the audio amplification device to the input thereof. 

25. A method for cancelling feedback in an audio amplification device 
comprising the steps of: 

applying an output of the audio amplification device to a first band 
limiting filter having a passband limited to a frequency band containing 
unstable frequencies. 

applying an output of the first band limiting filter to an adaptive 

filter; 

combining an output of the adaptive filter with an input of the 
audio amplification device. 

26. The method of claim 25 wherein the adaptive filter is implemented 
by a method comprising the steps of: 
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estimating power of a feedback cancellation error signal and an 
input signal to the adaptive filter; 

adapting filter coefficients in accordance with the estimated power; 
removing a DC offset from the adapted filter coefficients; 
bandpass filtering the adapted filter coefficients; 
applying the adapted filter coefficients to an adjustable filter. 

The method of claim 26 wherein the step of removing a DC offset 
less frequently than the step of adapting filter coefficients. 

28. The method of claim 26 wherein the step of bandpass filtering the 
adapted filter coefficients is performed only when an adapted filter coefficient 
exceeds a predetermined threshold. 

29. The method of claim 26 wherein the step of adapting filter 
coefficients comprises shifting, negation and addition without multiplication. 

30. The method of claim 25 wherein the adaptive filter is implemented 

by a method comprising the step of convolving the output of the bandpass filter 
with filter coefficients. 

3 1 . The method of claim 30 wherein the method of implementing the 
adaptive filter further comprises the step of applying a scaling gain to the filter 
coefficients. 

32. The method of claim 3 1 wherein the scaling gain is a power of 2. 


27. 

is performed 
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33. The method of claim 32 wherein the scaling gain is in the range of 2" 

3 to2 3 . 


34. The method of claim 25 wherein the audio amplification device 
comprises a hearing aid amplifier. 

35. The method of claim 34 further comprising the step of measuring 
adaptive filter coefficients with the hearing aid inserted in an ear of a wearer to 
identify unstable frequencies. 

36. The method of claim 25 wherein the first band limiting filter is a 
high pass filter. 

37. • The method of claim 36 wherein the passband of the high pass 
filter has a cutoff frequency approximately 200 Hz below a lowest of the unstable 
frequencies. 

38. The method of claim 25 further comprising the step of delaying the 

output of the audio amplification device prior to application to the first band 
limiting filter. 

39. The method of claim 38 wherein the output of the audio 
amplification device is delayed by an amount of time that is a function of a 
feedback path delay. 

40. The method of claim 25 further comprising the step of applying 
automatic gain control (AGC) to the input of the audio amplification device. 
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41 . The method of claim 40 wherein AGC is applied after the output 
of the adaptive filter is combined with the input of the audio amplification device. 

42. The method of claim 40 wherein AGC is applied to both the input 
of the audio amplification device and the output of the adaptive filter. 

43. The method of claim 25 further comprising the step of adjusting an 
adaptation speed of the adaptive filter. 

44. The method of claim 43 wherein the adaptation speed is adjusted as 
a function of a gain of the audio amplification device. 

45. • The method of claim 25 further comprising the step of selecting a 
length of the adaptive filter. 

46. The method of claim 45 wherein the length of the adaptive filter is 
selected as a function of sampling rate and: the frequency band containing unstable: 
frequencies so that duration of the adaptive filter's impulse response is long 
enough to cover an expected range of variation of group delay in the frequency 
band containing unstable frequencies in a feedback path from the output of the 
audio amplification device to the input thereof. 

47. The method of claim 25 further comprising the step of matching a 
frequency response of the first band limiting filter to a frequency response of a 
feedback path from the output of the audio amplification device to the input 
thereof. 
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18. A feedback canceller for an audio amplification device comprising: 

means for creating a first delay having an input coupled to an audio 
output of a hearing aid circuit and an output; 

a first band limiting filter having an input coupled to the output of 
the first delay means and an output; 

an adaptive filter having an input coupled to the output of the first 
band limiting filter and an output; 

means for creating a second delay having an input coupled to a 
conditioned output of a hearing aid microphone and an output; 

a first summing node having a non-inverting input coupled to the 
output of the second delay means, an inverting input coupled to the 
output of the adaptive filter and an output coupled to the input of the 
hearing aid processing module; 

a second band limiting filter having an input coupled to the input of 
the second delay means and an output; 

a second summing node having a non-inverting input coupled to the 
output of the second band limiting filter r an inverting input coupled to the 
output of the adaptive filter and an output; 

means for selecting a filter coefficient having a first input coupled 
to the output of the first band limiting filter, a second input coupled to the 
output of the second summing node and an output for supplying the filter 
coefficient to the adaptive filter; 

wherein the first and second band limiting filters have passbands 
limited to a frequency band containing unstable frequencies. 
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49. A feedback canceller circuit for an audio amplification device 
comprising: 

means for creating a delay having an input coupled to an audio 
output of a hearing aid circuit and an output; 

a first band limiting filter having an input coupled to the output of 
the delay means and an output; 

an adaptive filter having an input coupled to the output of the first 
band limiting filter and an output; 

a summing node having a non-inverting input coupled to a 
conditioned output of a hearing aid microphone, an inverting input coupled 
to the output of the adaptive filter and an output coupled to the input of 
the hearing aid processing module; 

' a second band limiting filter having an input coupled to the output 
of the summing node and an output; 

a third band limiting filter having an input coupled to the output of 
the first band limiting filter and an output; 

means for selecting a filter coefficient having a first input coupled 
to the output of the second band limiting filter and a second input coupled 
to the output of the third band limiting filter and an output for supplying 
the filter coefficient to the adaptive filter; 

wherein the first, second and third band limiting filters have 
passbands limited to a frequency band containing unstable frequencies. 

50. The device of claim 49 wherein the second and third band limiting 
filters have matching phase responses. 
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5 1 . The device of claim 49 wherein the second and third band limiting 
filters have substantially identical characteristics. 

52. The device of claim 1 further comprising means for storing a 
plurality of sets of filter coefficients for the first band limiting filter. 

53. The device of claim 52 wherein a set of filter coefficients is selected 
as a function of a gain of the audio amplification device. 

54. The method of claim 25 further comprising the step of storing a 
plurality of sets of filter coefficients for the first band limiting filter. 

55. The method of claim 54 further comprising the step of selecting a 
set of filter coefficients as a function of a gain of the audio amplification device. 


56. The method of claim 26 wherein the step of bandpass filtering the 

adapted filter coefficients applies a passband limited to a frequency band 
containing unstable frequencies. 
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